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Abstract 

The method for calculating the ground-state energy and the optical conductivity spectra is 
developed for a system of a finite number of interacting arbitrary-coupling polarons in a spherical 
quantum dot with a parabolic confinement potential. The path-integral formalism for identical 
particles is used in order to take into account the fermion statistics. Using a generalization of 
the Jensen-Feynman variational principle, the ground-state energy of a confined 7V-polaron system 
is analyzed as a function of N and of the electron-phonon coupling strength. The calculated 
optical conductivity spectra of the ./V-polaron system in a quantum dot manifest features related 
to ground-state transitions between states with different total spin. 

PACS numbers: 78.20.-e,71.38.+i,71.45.-d,71.45.Gm 
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I. INTRODUCTION 



Many-electron states in quantum dots are intensely investigated by various theoretical 
methods (see, e. g., Refs. [E 0, 0, HI ) • These works do not take into account the electron- 
phonon interaction, whereas it can contribute significantly to both the equilibrium and the 
non-equilibrium properties of quantum dots. In particular, the electron-phonon interaction 
plays a key role in the optical spectra of some quantum dots (see, e. g., Ref. p and 
references therein). To the best of our knowledge, the manv-polaron optical absorption in 
quantum dots has not yet been widely studied. In Refs. p, |7(, the ground state and the 
optical response of a fixed number of identical interactingpolarons are analyzed using the 
variational path-integral method for identical particles jaM 3| • As far as we investigate a 
system with a fixed number of identical particles, it should be described using the canonical 
ensemble As distinct from Ref. jfj, where only closed-shell systems are considered, the 
present approach is developed for both closed-shell and open-shell systems. 



II. THEORETICAL APPROACH 

We consider a system of N electrons with Coulomb repulsion, which interact with the lat- 
tice vibrations. We assume a spherical quantum dot with a parabolic confinement potential 
characterized by the frequency f2 . The partition function of the system is given by a path 
integral over electron and phonon coordinates (for details, see Ref. Q). After elimination 
of the phonon paths, the partition function of the electron-phonon system factorizes into a 
product of a free-phonon partition function with a partition function Z p ({iVo-} [where 
(3 = 1/ (ksT)] of interacting polarons, which is a path integral over the electron coordinates 
only: 



Px 

dx I D5c(r)e- S ^ T)] , (1) 



where S p [x (r)] contains the so-called influence phase of the phonons. It describes the 
phonon-induced retarded interaction between the electrons, including the retarded self- 
interaction of each electron with itself. The free energy of a system of interacting polarons 
F p {{N a } ,/3) is related to the partition function (JTJ) by the expression 

F p ({K} ,/3) = ~\nZ p ({N a } , 0) . (2) 

At present no method is known to calculate the non-gaussian path integral (H| analyti- 
cally. For distinguishable particles, the Jensen- Feynman variational principle [ll| provides 
a powerful approximation technique. It yields a lower bound to the partition function, and 
hence an upper bound to the free energy. 

It is a non-trivial problem, how a variational principle for the free energy should be 
formulated for a system of identical particles. However, it was shown that the path- 
integral approach to the many-body problem for a fixed number of N identical particles could 
be formulated in terms of a Feynman-Kac functional on a state space for N indistinguishable 
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particles, by imposing an ordering on the configuration space, and the introduction of a set 
of boundary conditions at the boundaries of this state space. The path integral (with the 
imaginary-time variable) for identical particles was shown to be positive within this state 
space. Thus a many-body extension of the Jensen- Fevnman inequality was found, which 
could be used for interacting identical particles (Ref. |8|, p. 4476). A more detailed analysis 
of this variational principle for both local and retarded interactions can be found in Ref. 
[Io| . It is required that the potentials are symmetric with respect to all permutations of 
the particle positions, and that both the exact propagator and the model propagator are 
antisymmetric (for fermions) with respect to permutations of any two electrons at any time. 
This means that both propagators have to be defined on the same configuration space. Under 
these requirements, the variational inequality for identical particles has the same form as 
that for the Jensen-Feynman variational principle: 

F p ^F + ±(S p -S ) So , (3) 

where Sq is a model action with the corresponding free energy Fq. The angular brackets 
mean a weighted ave-rage over the paths 

1 ^ (-1)^ 



((*)) So 7 UN \ fflE 



Z p ({N a },(3)^N 1/2 \N_ 1/2 \ 



X 



Fx 



rfx / Dx(t) (•)e~ 5o[s(T)] . (4) 



In the present work, we have chosen a model system consisting of N electrons with coordi- 
nates x = {x JiCr } and Nf fictitious particles with coordinates y = {yj} in a harmonic con- 
finement potential with elastic interparticle interactions as studied in Ref. 0. The model 
confinement frequencies for an electron and for a fictitious particle, the force constants and 
the mass of a fictitious particle are the variational parameters. 

In order to investigate the optical properties of the many-polaron system, we extend in 
Ref. 0] the memory-function formalism developed in Ref. jl2| to the case of interacting 
polarons in a quantum dot. Within this technique, the optical conductivity for a system of 
interacting polarons in a parabolic confinement potential is given in terms of the memory 
function x (^)j 

Rea(u) = 5 5-, (5) 

m[^-^-RexM] 2 + [Im X H] 2 ' 

where m is the electron band mass, and x ( w ) is 



2.2 CV 



X ( U )=Y f (e**-l) 



o 



xIm[T: LO (t)( Pq (t)p_ q (0)) M ] dt. (6) 

Here, Vq is the amplitude of the electron-phonon interaction [l2| , (t) = 
cos [uj (t — iH/3/2)] I sinh (/3hu>/2) is the phonon Green's function, and (p q (t) p_ q (0)) M is 
the density-density correlation function calculated using the model Lagrangian of electrons 
harmonically interacting with fictitious particles. It is noteworthy that the optical conduc- 
tivity © differs from that for a translationally invariant polaron system both by the explicit 
form of x (w) and by the presence of the term Qq in the denominator. 
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III. RESULTS AND DISCUSSION 



Further on, we present numerical results for the many-polaron system in a quantum dot 
at T = 0. The total spin of interacting polarons in its ground state, as well as the total 
spin of a confined few-electron system, depends on the number of electrons. The analysis 
of this dependence reveals a shell structure of quantum dots. In Fig. Q the total spin S 
of a finite number of interacting polarons in a spherical quantum dot for different values of 
the confinement energy ftflo, of the electron-phonon coupling constant a and of the ratio of 
the high-frequency and the static dielectric constants r\ = Eoo/eq is plotted as a function of 
the number of electrons in a quantum dot. As distinct from few-electron systems without 
the electron-phonon interaction, three types of spin polarization are possible for the ground 
state, which can be distinguished from each other using, e. g., capacity measurements. 




FIG. 1: Total spin of the system of interacting polarons in a quantum dot as a function of the 
number of electrons for the confinement energy H£Iq = 0.5H* (a) and for htlo = 0.1H* (b). The 
confinement energy is measured in effective Hartrees H* = \m/ (jnos^)] Hartree, where mo is the 
electron mass in the vacuum. 



(i) Except the strong-coupling regime and the low-density case, for closed-shell systems 
5 = 0, while for open-shell systems S takes a maximal value for a given shell filling, in 
accordance with the Hund's rule. 
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(ii) When weakening confinement for a fixed number of electrons, the electron density 
lowers. Hence, at sufficiently small values of Qq, a spin-polarized state for a system of 
interacting polarons in a quantum dot becomes energetically more favorable than a state 
satisfying the Hund's rule. 

(iii) In the strong-coupling regime (a> 1 and rj <C 1), the total spin of an open-shell 
system for the ground state can take a minimal possible value. This trend to minimize the 
total spin is a consequence of the electron-phonon interaction, presumably due to the fact 
that the phonon-mediated electron-electron attraction overcomes the Coulomb repulsion. 

The optimal values of the variational parameters are used as input for the calculation of 
the optical-conductivity spectrum of the system (JSJ. The ground-state transitions between 
the aforesaid states with different values of the total spin, which occur when changing the 
number of electrons, are manifested in the optical absorption spectra. In particular, the 
shell structure for a system of interacting polarons in a quantum dot is clearly pronounced 
when analyzing the first frequency moment of the optical conductivity, 
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FIG. 2: The function (N) (a) and the addition energy A (iV) (6) for systems of interacting 
polarons in a quantum dot with a = 3, rj = 0.3 and f^o = 0.5cjlo (^o ~ 0.01361 H*). 

Figure El represents the function 

B(N)= (u)\ N+1 -2{u)\ N +{u;)\ N _ 1 (8) 
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along with the addition energy A (JV) needed to put an extra electron into the quantum dot 
with N electrons, 

A (iV) = E° (N + 1) - 2E° (TV) + E° (N — 1) , (9) 

where _E° (AT) is the ground-state energy. 

Both functions 6 (N) and A (N) exhibit features at A" = 4 and at N = 7, where the 
aforesaid transition occurs. Namely, at A r = 4, © (N) and A (N) have minima, whereas 
at AT = 7, O (N) has a minimum, and A (N) has a kink. Pronounced peaks appear in 
these functions at the "magic numbers" A^ = 10 and A" = 20, which correspond to closed- 
shell systems. The transition between the spin-polarized ground state and the ground state 
obeying the Hund's rule thus should be observable using optical measurements. 



IV. CONCLUSIONS 

The path-integral treatment of the quantum statistics of indistinguishable particles al- 
lowed us to find an upper bound to the ground-state energy of a fixed number of polarons 
in a parabolic confinement potential for arbitrary electron-phonon coupling strength. 

The transitions between states with different values of the total spin manifest themselves 
through discontinuous changes in the optical conductivity spectra and of the addition energy 
as a function of the number of electrons. Thus, the analysis of the optical-conductivity 
spectra provides a tool for examining the shell structure of a system of interacting polarons 
through experimental measurements of the optical absorption in quantum dots. 
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